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Abstract 

The value distribution of derivatives of characteristic polynomials of matrices from 
SO{N) is calculated at the point 1, the symmetry point on the unit circle of the eigen- 
values of these matrices. We consider subsets of matrices from SO{N) that are con- 
strained to have n eigenvalues equal to 1, and investigate the first non-zero derivative 
of the characteristic polynomial at that point. The connection between the values of 
random matrix characteristic polynomials and values of L- functions in families has been 
well-established. The motivation for this work is the expectation that through this con- 
nection with L-functions derived from families of elliptic curves, and using the Birch 
and Swinnerton-Dyer conjecture to relate values of the L-functions to the rank of elliptic 
curves, random matrix theory will be useful in probing important questions concerning 
these ranks. 

1 Introduction 

1.1 Random matrix theory and number theory 

The connection between random matrix theory and number theory began with the work 
of Montgomery when he conjectured that the distribution of the complex zeros of the 
Riemann zeta function follows the same statistics as the eigenvalues of a random matrix 
chosen from U{N) generated uniformly with respect to Haar measure. This conjecture is 
supported by numerical evidence [H^ ^-Iso by further work |17[ I34[ [21 [H] suggesting that 
the same conjecture is true for more general L-functions. For all these L-functions there 
is a Generalized Riemann Hypothesis that the non-trivial zeros lie on a vertical line in the 
complex plane. The conjectures mentioned above concern the statistics of the zeros high 
on this critical line. 

The philosophy of Katz and Sarnak j2()| I21j extended the connection with random 
matrix theory by proposing that rather than averaging over many zeros of a given L- 
function, if the zeros near to the point where the critical line crosses the real axis are 
averaged over a family of naturally connected L-functions, then they will be found to follow 
the statistics of the eigenvalues of one of the three classical compact groups of random 
matrices: U{N), 0{N) or USp{2N), where again the statistics are computed with respect 



to the probability measure given by Haar measure. There is numerical evidence for this 
conjecture as well |38j . and strong support is given to it by the rigorous work of Katz and 
Sarnak in the case of function field zeta functions. 

For a review of applications of random matrix theory to questions in number theory 
see, for example, jl] or 

There has been a series of papers, starting with [22] and continuing with [HI I19| 
I18| I22| IHl Cj , examining how random matrix theory can be used to predict the distribution 
of values of the Riemann zeta function and other L-functions, either averaged over an 
interval high on the critical line, or over a family at the critical point where the critical line 
crosses the real axis. For large values of the natural asymptotic parameter, for example 
the variable ordering the L-functions within the family, the moments of L-functions are 
conjectured to split into a product of an arithmetic contribution, determined by the family 
being averaged over, and a component derived from a random matrix calculation - the 
corresponding moment of the characteristic polynomial of the matrices in one of the three 
matrix groups U (N), 0{N) or USp{2N). The asymptotic parameter on the random matrix 
side is the dimension of the matrix and a natural equivalence can be made between the 
two. 

As will be reviewed in the following section, this conjecture for the leading order 
behaviour of moments of L-functions has been used by Conrey, Keating, Rubinstein and 
Snaith [Hj to predict the frequency, within a family of elliptic curves, of an L-function taking 
the value zero at the critical point. This is done by using random matrix theory to predict 
the value distribution of the L-values at the critical point and then the discretization of 
these values (HHlESlllSl to calculate a probability of the L-value being zero. See also David, 
Fearnley and Kisilevsky J2j for a similar use of random matrix theory to predict frequency 
of vanishing at the critical point amongst families of elliptic curve L-functions twisted by 
cubic characters. 

In Section 12.11 we calculate the value distribution of the first non-zero derivative of 
the characteristic polynomial at the point one (corresponding to the value at the critical 
point in the analogy with L-functions) when a given number of eigenvalues are conditioned 
to lie at the point one. We find, at equation 1)2. 2U() . that moment of the nth derivative 
grows like 

M{n,M,s) := (|A[;^(l)r) ~ (n!)^(27r)^/^2-^'/2-s(n.-i) + 1/2) ^ .2/2+,(„-i/2)^ 

Ct{ti -\- 1/2 -|- sj 

(1.1) 

where Aij{e'^^) is the characteristic polynomial of U, G{z) is the Barnes double gamma 
function (see l2.l4|) and the angle brackets denote an average over the set of matrices from 
SO{N) with n eigenvalues lying at the point one (A'" = n + 2M). We also show that the 
probability that |A[]'^(1)| < X over those U G SO{N) with n eigenvalues at I (again, 
N = n + 2M) is, for small X, given by (see ()2.12() and the sentence following) 

2n+l 

^^X^f{n,M), (1.2) 

for the function / given at ()2.13|) . 

If a discretization of the values of derivatives of L-functions at the critical point 
were known, this would provide a way to predict the frequency of zeros of various orders 
occurring at this point. This has been investigated in the case of a single zero at the critical 
point JJ. In this situation, n = 1, the results of |llj give support to the validity of the 
model presented here. 
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The result presented in Section on the moments of the first derivative of charac- 
teristic polynomials from S0{2N + 1) (equation II. ll with n = 1) has already been applied by 
Delaunay in in order to predict the moments of the orders of Tate-Shafarevich groups 
and regulators of elliptic curves with odd rank belonging to a family of quadratic twists. 

In Section 12.21 we briefly discuss the eigenvalue statistics of random matrices with 
a certain number of eigenvalues conditioned to lie at one. 

We note that while we believe that the model presented in Section [2.11 should apply 
to the L- functions selected from families of quadratic twists of elliptic curves for the property 
that they have high order zeros at the critical point, there has been some very interesting 
theoretical work computing the one- and two-level densities for parametric families of elliptic 
curves that implies that the zeros of the associated L-functions follow a different model (see 
j28j and j4Uj ) . This does not seem to be a contradiction, as the zero statistics are examined 
over collections of L-functions selected in very different ways, but it certainly makes the 
question of the zero statistics of L-functions with high-order zeros at the critical point very 
intriguing. See Section |21 for further discussion on this issue. 



1.2 Random matrix theory and elliptic curves 

We review here the results of jHI which apply random matrix theory to predicting the 
frequency of vanishing at the critical point of the L-functions in the family of elliptic curves 
described below; or equivalently, assuming the Birch and Swinnerton-Dyer conjecture, the 
frequency of rank 2 or higher curves occurring in the family of elliptic curves. The motivation 
for the random matrix calculations presented in Section 12.11 is the expectation that they 
may be used similarly to examine curves of higher rank. 

We consider an L-function (defined by the Dirichlet series and Euler product below 
when Re s > 3/2) 

oo 

^^(') = E 5 = U^piW) = n (1 - '^pp-T'Ui^ - ^pp" +p'-'T' ' (1-3) 

n=l p p\A pfA 

that is associated to an elliptic curve E over Q 

E : y"^ = x^ + Ax + B. (1.4) 

The coefficients Cp, for prime p, are determined by = p -|- 1 — ^E{¥p), where i^E{¥p) 
counts the number of pairs x,y, with < x,y < p — 1, such that = -|- Ax + B ( 
mod p), plus one for the point at infinity. A is the discriminant of the cubic -|- Ax + B. 
For an extremely clear introduction to elliptic curves in the context discussed here, see the 
review paper by Rubin and Silverberg [5^ . 

A family of elliptic curves is formed by 

Ed : dy"^ = x^ + Ax + B (1.5) 

for integers d that are fundamental discriminants. (A fundamental discriminant is an integer 
other than 1 that is not divisible by the square of any prime other than two and that satisfies 
d = 1 mod 4 or (i = 8, 12 mod 16.) The corresponding family of L-functions, ordered by 
\d\, are 

M^,x.) = f;^, (1.6) 

n=l 
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where the characters Xdi'iT') = (^) are Kronecker's extension of the Legendre symbol, defined 
for prime p as 

+1 a p \ d and = d{modp) is soluble 
ifp\d . (1.7) 

— 1 if p \ d and = d{modp) is not soluble 

The family of curves is called the family of quadratic twists of E because the characters 
Xdin-) are real, quadratic Dirichlet characters. 

Evidence (see for example jSB]) points to the zeros near the critical point (where the 
critical line crosses the real axis) of such a family of L-functions having statistics like the 
eigenvalues near 1 of matrices from 0{N) with Haar measure. To explain this correspon- 
dence more clearly, about half of the L-functions in this family are expected to have an even 
functional equation that relates LE{s,Xd) to Le{2 — s,Xd), meaning that the zeros on the 
critical line (Res = 1) appear in complex conjugate pairs and if there is a zero at the critical 
point then it must be of even order, while the other half have an odd functional equation, 
forcing a zero of order 1,3,5,... at the critical point and all other zeros appearing in complex 
conjugate pairs. Matrices in 0{N) also have eigenvalues that appear in complex conjugate 
pairs, with the possible exception of unpaired zeros at 1 and -1. The low- lying zeros of 
elliptic curve L-functions from the family above that have even functional equation show 
the statistics that eigenvalues of S0{2N) display near the point one. These eigenvalues 
occur in complex conjugate pairs and an eigenvalue at one must have even multiplicity. In 
contrast, the low- lying zeros of the L-functions above that have an odd functional equation 
display the same statistics as SO{2N + 1) eigenvalues near one, since in this case there 
is always an eigenvalue at one and it has to have odd multiplicity. (Note that we could 
just as well have chosen 0~{2N) matrices (orthogonal with determinant -1) to model the 
L-functions with odd functional equation, since the statistics of these eigenvalues near to 
one are the same as those of SO{2N + 1), see for example |20j.^ 

Both numerical and analytical evidence has been given already |H] that random 
matrix theory can be used to conjecture the frequency of L-functions with even functional 
equation vanishing at the critical point in a family such as that described above. This is 
particularly important because of the Birch and Swinnerton-Dyer conjecture which suggests 
that the order of the zero of an elliptic curve L-function at the critical point is the same as 
the rank of the elliptic curve itself. The rank is a non-negative integer r that characterizes 
the number of rational points on the elliptic curve. The set of pairs x, y G Q that satisfy 
equation (|1.4|) . plus one point at infinity, form a commutative group (with a standard law 
of addition defined on the curve) that is isomorphic to Z''' ^ E{Q)tors- E{Q)tors is the 
subgroup of elements of finite order and r is the rank of the curve. The frequency of various 
ranks in a family such as the one defined at H1.5|) . or even whether the rank is bounded in 
such a family of elliptic curves, are important unanswered questions. 

In inj the first step is taken towards using random matrix theory to generate con- 
jectural answers to these questions by addressing the frequency of L-functions from elliptic 
curve families having any zero at all at the critical point. The argument used will be 
sketched below to illustrate the role played by random matrix theory. The purpose of this 
paper is to perform the random matrix calculations needed to predict the frequency of a 
zero of a given order at the critical point. Applying the random matrix results from Section 
12.11 to the order of vanishing of L-functions involves the knowledge of complicated arith- 
metic quantities and much more work is needed in this area. Some preliminary numerical 
investigation has been done for the case of third order vanishing. 




4 



We now review the results in [HI- We define the family of L- functions discussed 
above with even functional equation: 

Te+ = {LE{siXd) having even functional equation}. (1-8) 

The zeros of these L-functions near the critical point, s = 1, show the same statistics as the 
eigenvalues of matrices from SO{2N), and so the L- values at the critical point are modelled 
by the characteristic polynomial 

TV 

Au(e'') = n (l - ^'^'""'0 " ' (1-9) 

n=l 

evaluated at the point ^ = 0. Here e^*^^ , • • • , e^*^^ are the eigenvalues of the matrix 
U G S0{2N). 

The moments of Ac/(1) = n^=i 1^ ~ | easily calculated using Weyl's ex- 
pression for Haar measure on the conjugacy classes of S0{2N) and a form of Selberg's 
integral to be [221 



/ AuiiydUHaar 
JS0(2N) 

_ ^2Ns fr r(iv + j-i)r(g + j-i/2) 

" /jr(j-i/2)r(s+j + iv-i) ^'-^'^ 

= Mo{N,s). 

The L-function moments are then conjectured to have the form pi I22j 



Me(r,s) = — Yl LEihXdf 

\d\<T 

- as{E)Mo{N,s) (1.11) 

for large T. Here = log T (from equating the density of zeros near the critical point with 
the density of the matrix eigenvalues), the sum is over fundamental discriminants d, T* is 
the number of terms in the sum and as{E) is an Euler product that contains arithmetic 
information specific to the elliptic curve E and the family of L-functions being averaged 
over. In practice, it is often a subset of J-e+ that is summed over. If, for example, we select 
those L-functions LE{s,Xd) in with d > and further restricted by a condition on d 
mod Q, if Q is odd, and on d mod 4(5, if Q is even, then the arithmetic factor would be 



p\Q 

See and [HI, Section 4.4, for more examples. 



5 



Next we consider the distribution of the values of the characteristic polynomials of 
SO{2N) matrices at the point 1. If Po{N,x)dx is the probability that the characteristic 
polynomial of a matrix chosen from SO{2N) with Haar measure has a value between x and 
X + dx, then 

Po{N,x) = ^ [ Mo{N,s)x~'ds 
2mx 

~ x~^/'^h{N) (1.13) 

for X — > 0"*", since for small x the behaviour is dominated by the pole of Mo{N,s) at 
s = —1/2. Here (c) denotes a path of integration along the vertical line from c — ioo to 
c + ioo. 

For large iV, h{N) ~ 2~'^/^G{1/2)tt~'^/'^N^/^ {G is the Barnes double gamma func- 
tion, defined as [H l37j: 

oo 

G(l + z) = (27r)^/2e-[(i+^)^'+^]/2 J] [(1 + z/?i)"e-^+^'/(2n)j ^ 

n=l 

where 7 is Euler's constant. See also (|2.14|) for more properties of this function.) Since the 
probability that an element of S0{2N) has a characteristic polynomial whose value at 1 is 
X or smaller is Po{N, x)dx, we find that the the behaviour of this probability for small 
X and large N is 

rX 

lim lijn / Po(N,x)dx) =2^/^G(l/2)Tr~^/'^. (1.15) 

AT^oo X-*0+ ^ Jo 

We see from equation (jLlip that for large d, moments of L- functions are conjectured 
to be just as{E) (the prime product) times the random matrix moment Mo{N, s). If this is 
true, then we define Pe(T, x)dx as the probability, amongst members of Te+ j that Le(1, Xd), 
for \d\ around e^, will take a value between x and x + dx, giving 

Pe{T,x) = ^ [ ME{T,s)x-'ds, (1.16) 
2mx 

and an approximation for this probability for small x should be 



Pe{T,x) ~ a_,^2{E)x-^/^h{N). (1.17) 

Here equating densities of zeros means N ~ log T. 

So, roughly, random matrix theory predicts that the probability that an L-function 
LE{s,Xd) G ^E+j with \d\ close to t, has a value at s = 1 that is X or smaller is 

~ 2a_i/2^i/22-^/8G(l/2)^-i/4 log3/8 1. (1.18) 

But these L-functions are constrained to take only certain discretized values at the 
critical point s = 1. The L- values have the form [38l l35l [25] : 

LE{l,Xd) = ^E^^^, (1.19) 
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where the C£;(|d|) are integers, the Fourier coefficients of a half-integral weight form. 
So, if 

LE{l,Xd)<^ (1.20) 
\/\d\ 



then 

LE{l,Xd)=^ (1.21) 

and the first thought would be to take X = KEl\/t in (|1.18|) and then to integrate t from 
to T . However, there is arithmetical information encoded in the C£;(|d|) (one is not always 
the smallest allowed non-zero value) and so to avoid this problem, the conjecture stated in 
[H] is restricted to prime discriminants. 

Conjecture 1.1 (Conrey, Keating, Rubinstein, Snaith): 

Let E he an elliptic curve defined over Q. Then there is a constant ce > such that 



^ 1 CET''/\logTy 



5/8 



iB{l,Xp)=0 



(The conjecture was originally stated in [Hj with ce > 0, but in fact numerics have revealed 
that the constant can be zero ^01 for certain families. An explanation of such a case with 
ce = was given by Delaunay Jl].) 

With the Birch and Swinnerton-Dyer conjecture, this suggests that out of all the 
elliptic curves associated with L-functions in J-'e+ with prime discriminant p <T (there are 
of order T/ log T of them), a number of order r^/'^(log T)~^/^ of them should have rank two 
or greater. The T^/^ has been predicted previously by Sarnak using different arguments, 
but random matrix theory adds more detailed information in the form of the power on the 
logarithm. For numerical evidence supporting the conjecture, see jlj and fTO] . 

In the next section we will calculate the distribution of values equivalent to H1.13() 
that is required to probe questions of elliptic curves of a given rank occurring in families of 
quadratic twists (equation 1)1. SI) ). 



2 Random matrix calculations 

In this section we will calculate the probability of the nth derivative of a characteristic 
polynomial ()2.1)) of a random SO{N) matrix taking a value less then X at 6 = when all 
lower derivatives are constrained to be zero at this point. That is, we perform the random 
matrix average over the (measure zero) subset of SO{N) of matrices with n eigenvalues at 
the point 1 on the unit circle, and 2M other eigenvalues occurring in complex conjugate 
pairs {N = n + 2M). 

We will also illustrate the n-level density of the 2M eigenvalues not constrained to 

lie at 1. 

2.1 Value distribution of higher derivatives 

For a matrix U G SO{N) with n eigenvalues equal to one (A^ = n + 2M), the characteristic 
polynomial looks like 

M 

Auie'^) = (1 - e"^^)" JJ(1 - e*(^^-^))(l - e^^"^^"^)). (2.1) 
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We consider the nth derivative 

M 



-(1 - e-^r - e^^^"")(l - e-^^-°)|„^, 



da 



4— (1 -e-")" 



:0 



and use this to define 

M M 

a|;)(1) = n! ]^(1 - e*^0(l " e"^^0 = ^!2^'^ " ^"^^i)' (2-2) 

The matrices with n eigenvalues equal to one form a set of measure zero within 
SO{N) with Haar measure, for n > 1. For even = 2P, Haar measure on the conjugacy 
classes of S0{2P) is given by 

C JJ {cos9k-cos9jfdei---dep, (2.3) 

i<j<k<P 

with C a normalization constant, and we define the probability of having 2q eigenvalues at 
1 and the rest within infinitesimal neighbourhoods of e^*^^ , ■ ■ ■ , e^*^^~« as 

(lo ■■■Jo ni<j<fc<p(cos 9k - cos 9jfd9p-q+i ■ ■ ■ d9p\ d9i--- d9p-q 

liin 

!o ■■■ lo /o° • • • lo ni<j<fc<p(cos 9k - cos 9jYd9p-q+i ■ ■ ■ d9pd9i ■ ■ ■ d9p_g 

= C{P,q)\\_{l-cos9jf'i {cos9j -cos9kfd9i---d9p_g. (2.4) 

j = l l<j<k<P-q 

Here C{P,q) is the normalization constant. 

The same procedure for N odd, starting from the measure on S0{2P + 1) 
p 

C" Ylil- cos 9 j) Yl {cos9k-cos9jfd9i---d9p, (2.5) 

j=l l<j<k<P 

with normalization constant C", leads to a general expression for the measure on the matri- 
ces in SO{N) (where now N can be even or odd) with n eigenvalues equal to 1 {N = n+2M): 

M 

C(M,n)J]( l-cos6'j)" JJ (cos6'j- - cos6'fc)^. (2.6) 

j=l l<j<k<M 

Thus the quantity to calculate is 
M{n,M,s) := 

M 



c{M,n) ... iA[;)(i)rn( 

1 - cos 6*^)" (cos 9j — cos 9kfd9i • • • d9M 

j=l l<j<k<M 

= C{M,n) f ... f (n!)^2*^"]^(l -cosej)"+" {cos9j - cos9kfd9i . . . d9M 

•^0 j=l l<j<k<M 

7 •••/ n n ^ u/2 n i^j-^k?dx,...dxM.{2.7) 



C7(M, n)(n!)"2 
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This can be evaluated using a form of Selberg's integral (for details see j26jl: 



-1 rl K 



I, '"I, n \i^j-^i)\''tlii-xjr-Hi+xjf-'dxj 



^<j<l<K 3=1 
^ ^yKiK-l)+K{a+/3-l) JT r(l + 7 + jj)r{a + jj)r{(3 + j-f) 

j}^ r(i + 7)r(a + /3 + 7(K + j-i))' ^ • 

if Rea > 0, Re/3 > and Re7 > - min (;^, frr, 

We have 7 = 1, a = n + 1/2 + s and /? = 1/2, so the integral in equation (|2.7p is 

M{n,M,s) = C{M,n){n\y2^^''+^'^^'-^^+^^ 

^ ^ r(j + i)r(n-i/2 + g + j)r(i-i/2) 

r(n + M + s + j-l) ■ 
But the normalization constant can also be evaluated by Selberg's integral to be C{M, n) = 

JLir(n-l/2+j)r(n + s + M+j-l) ^ ^ 

Consider n fixed and finite and M large. If the probability that |a[]^^(1)| takes 
a value between x and x + dx is given by -P(n, M, x)dx then from a standard result in 
probability (with (c) denoting a path of integration along the vertical line from c — ioo to 

c + ioo) 

P(n,M,x) = ^— / x''M(n,M,s)ds (2.11) 
27rzx ./(c) 



(c) ^ f}j{n-l/2+mn + s + M + j-ir- 



2'Kix 



We are interested in the behaviour at small and this is dominated by the nearest 
pole to zero of the integrand: the pole at s = — (n + 1/2) of T{n — 1/2 + s + j). Thus for 
small X 

P{n, M, x) ~ M). (2.12) 

The probability that lAf^^l)! < ^ 

over those U G SO{N) with n eigenvalues at 1 is 

therefore ~ 2n+i -^ ^ ^'^^ small X and the dependency on n and M is determined 

by 

-2n-l w 2n+l s 1 il r(o-)r(M + n + 7 - 1) 

We can discover the behaviour of /(n, M) for large M by introducing the Barnes 
G-function 

00 

G{1 + z) = (2^)-/2e-[(i+7)-'+-]/2 Yl [(1 + z/n)"e-^+^'/(2n)J ^ (2.14) 

n=l 
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5 10 15 20 



P(X) P(X) 
O.lr O.lr 




Figure 1: Distribution of values {P{n, M, x) from 1)2.11(1 ) at the point one of the nth deriva- 
tive of the characteristic polynomial of the subset of matrices from SO{n + 2M) that are 
conditioned to have n eigenvalues at one. Figures c) and e) show in more detail the be- 
haviour at the origin of figures b) and d) (see equation (|2.12j) ). 
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which has zeros at the negative integers, —re, with multiphcity n (n = 1,2,3...). Other 
properties useful to us are that 



G(l) = 1, (2.15) 
G{z + l) = r{z)Giz), 

and furthermore, for large \z\ 

log G{z + 1) ~ ^2(1 log z-l) + ^z log(27r) - ^ log z + C'(-l) + 0(i). (2.16) 
Thus 

So we can write 

f(n M) - (rJ)"^2-^^(^^ G{M)G{2M + re)G(re + l/2)G(Af - 1/2) 
f{n,M)-[n.) 2 G(M + re)G(re + M + 1/2)G(2M - 1/2) ^^'^^^ 

and expanding the G-functions for large M gives 

— 2n— 1 n 3 3ri 7 n 1 

f{n,M) ~ (n!)^~G(re + l/2)M~~+2 + 82~"~~87r~2-4. (2.19) 

Note that we can also use ()2.16|) to revisit M.{n, M, s) (defined at l2.7|) and examine 
that asymptotically for large M. This gives us the large M behaviour of the moments at the 
point one of the reth derivative of characteristic polynomials for which all lower derivatives 
vanish. We have 

M(„ M = (^^Y-JMs G{n + l/2 + s + M) G(n + 1/2) G(re + 2M) G{n + s + M) 
^ G{n + 1/2 + s) G{n + l/2 + M) G{n + M) G{n + s + 2M) 

~ (n!)«(27r)^/22-^'/2-s(n-i)_^(!l±V^^sV2+s(n-i/2)_ ^2.20) 

G{n + 1/2 + s) 

2.2 Eigenvalue statistics 

To get an idea of what the eigenvalue statistics of the remaining eigenvalues look like when 
we condition re eigenvalues to lie at one in the manner described in (|2.4|) . we briefly review 
the m-level densities of the unconditioned eigenvalues. These calculations are not new; 
these same matrix ensembles were studied by Nagao and Wadati [30] and also Duenez 
jl5j and are discussed in the book by Forrester JH]- Concurrently with the present work, 
this matrix model has also been discussed in connection with the zero statistics of elliptic 
curve L-functions by Miller (with appendix by Duenez). However, we include a brief 
discussion of the statistics here to illustrate the effect that forced eigenvalues at one have 
on the remaining eigenvalues. 
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The terminology for these statistics varies (for example, Mehta [^H] refers to them 
as m-point correlation functions) so we will write (with C the normalization constant) 



cf.-.r Yi /(%,...,%„) (1-COS0 

° ° i<ii<-<jm<A' 

(cos 9j — cos Oki^dOi • • ■ cIOm 



X 

1 r r 

ml Jo Jo 



f{6i, . . . , 9m) 

m 

X det{K^lf^'^'-^'^\d,,ek)),,k=i,...,m 11(1 - cosd.Tddi ■■■d9m 

• •• r f{ei,...,em) (2.21) 

Jo 

X det(2" sin" | sin" % i^Sr'^''"'^'Vi, ^fc))i,fc=i,...,m^^ei • • • 



and refer to det(2"sin" ^ sin" ^ i^T^^ ^^^^ (%, ^fc))j,fc=i,...,m as the m-level density (see 
[H] for further discussion of these densities). By standard techniques (see for example j2S]; 
chapter 5, [30] i chapter 4) 

{fj,"k))j,k=l,...,m 

M 

= |;(/if V'^^""'^''"'^'^ (cos ^,) P,('^-^/^'-^/^)(cos^fc) 

j=0 

2-n+i r(M + 2)r(M + n + l) 



2M + n + 1 r(M + n + l/2)r(M + 1/2) 



(2.22) 



/p(n-l/2,~l/2). . X p(n-l/2,-l/2). . x p(n-l/2-l/2). . x p(n-l/2-l/2). . 
^ / (cos^jj l^jy^ (cos dk) - t^M [COSUj) (COS f^k) 

I COS 9j — cos 9k 

is defined in terms of the Jacobi polynomials Pj°''^\cos9) (described in detail in [HS]) with 
the orthogonality condition 

Pj"'^^(cos0) P^"''^^(cos0)(l - cos0)"d0 = Sjkhl^'f^\ (2.23) 







where 

hp^' = — - y ' - ' ' ^ ' ^( ■ (2.24) 



^{a,(3) ^ 2^+^+^ r(j + g + i)r(j + P + 

2j + a + /3 + 1 r(j + l)r(j + a + /3 + 1) 



Scaling the eigenvalues and using asymptotic formulae for the Jacobi polynomials 
(see [SSj, page 197) we find that in the large- matrix limit (where Ja{x) are Bessel functions 
of the first kind) 

^(n-i/2,-i/2)(^^^) ji^n^^j4-V2,-i/2)(^^^)2"sin"(^)sin"(^) (2.25) 



/ /,Nl/2/ ,^l/2^'^n~3/2(^7^)J„_l/2('^7r) - (^J„_3/2((^7r)J„_i/2(6l7r) 

{tt9) (vrc/.) 
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The scaled 1-level density is therefore 

^(n-lA-l/2)(^^^) = y^(J^3/2(^^) + ^n-l/2(^^)-^^"/n-l/2(^^)^„-3/2(evr)), (2.26) 

and can be seen to behave hke a constant times for small 9, where we remember that 
n is the number of degenerate eigenvalues. This increasing repulsion of the first eigenvalue 
not located at 1 can be seen in Figure [21 




Figure 2: Scaled 1-level densities, from equation (|2.26j) for various values of n, the degener- 
acy of the eigenvalue at one. 
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3 Discussion 



We have shown that the probabihty that |a[^^(1)| < X over those U G SO{N) with n 
eigenvalues at 1 (N = n + 2M) is, for small X, given by 

2n+l 

^^X^f{n,M) (3.1) 
and the mean value of the nth derivative grows like 

(jr\Ti -\- 1/2 -|- S) 

(3.2) 

where the angle brackets denote an average over the set of matrices from SO{N) with n 
eigenvalues lying at the point one (A^ = n + 2M). 

The result (|3.2|) , see also (j2.2()|) , was applied by Delaunay in ^3] to predict moments 
of the orders of Tate-Shafarevich groups and the regulators of elliptic curves belonging to 
a family of quadratic twists, see (|1.5|) . The Birch and Swinnerton-Dyer conjecture provides 
a formula for the first non-zero derivative of an L-function in terms of various quantities 
related to the associated elliptic curve. One of these quantities is the order of the Tate- 
Shafarevich group, and another, in the case of the first derivative of L-functions with odd 
functional equation, is the regulator. For families of L-functions with even functional equa- 
tion the result ()2.2U|) is used with n = by Delaunay to predict the asymptotic form of 
moments of the order of the Tate-Shafarevich group for the associated family of elliptic 
curves, and for families with odd functional equation the case n = 1 was used to conjecture 
the form of moments of the regulator. 

To predict, for example, the number of L-functions in T^- (the family defined like 
J-E+ in (|1.H() but with odd functional equation) that have a zero of order at least three at 
the point 1 we need the result (|3.1j) for n = 1 plus some information analogous to (|1.2()|) 
giving the smallest non-zero value that can be taken by the derivative of an L-function 



from this family. A formula like (|1.20j) is not known at the moment for derivatives of 
L-functions, but numerics suggest that there is a gap between zero and the smallest 
non-zero value of L', which allows random matrix theory to be applied to the question of 
the number of L-functions with a zero of order at least three at the point 1. Further work 
in this area is ongoing, but it is also worth noting that the numerics in support the 
random matrix model presented here, as the numerical cumulative distribution of values of 
L' for L-functions in J^^- shows the expected x^/^ behaviour near x = 0, agreeing with the 
exponent 3/2 in (|3.H) when n = 1. 

This leads us to believe that the measure in 1)2. 6() should correctly model the zeros 
near the point 1 of L-functions with at least an nth order zero at 1 selected from the 
family J^^- or (depending on whether n is odd or even, respectively), but we note 
again that this is not the only model for L-functions with n zeros at 1. The measure 
presented here is what Miller [27j calls the Interaction Model, and he contrasts this with 
the Independent Model where the n eigenvalues at 1 have no effect on the statistics of 
the remaining eigenvalues. This gives, for example for matrices from S0{2M + r) with r 
eigenvalues at 1, the set of matrices 

I J. X J. 



) :9G50(2M)| 
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with a r X r identity block in the upper left corner and with the joint probability density 
on the remaining eigenvalues being Haar measure on SO{2M), that is 

cx JJ {cos 9k — cos 6 j)'^ JJ dOj. 
l<j<k<M i<j<M 

Miller finds numerical evidence j27j . and for theoretical results see HSj and j4Uj . that the 
Independent Model agrees with the zero statistics of L-functions of a parametric family of 
elliptic curves constructed so that each member of the family has rank r. Through the 
Birch and Swinnerton-Dyer conjecture, the associated L-functions are expected to have a 
zero of order r at the point 1, and Miller finds evidence for this high-order zero, as well 
as for the fact that its presence does not affect the statistics of the other zeros in the way 
that the model presented here would (for example as illustrated in Figure 2), but rather the 
zeros follow Miller's Independent Model. This is not a contradiction to the proposal that 
the model presented here (the same as Miller's Interaction Model) is correct for the family 
of quadratic twists of elliptic curve L-functions. It will be the subject of further work to 
produce convincing numerics for the zero statistics of the elliptic curve L-function families 
discussed in this paper, although this is a difficult task due to the small number of curves 
in these families with rank greater than or equal to two. 
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